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ABSTRACT 

To compare photometric properties of galaxies at different redshifts, the fluxes need 
to be corrected for the changes of effective rest-frame wavelengths of filter bandpasses, 
called if -corrections. Usual approaches to compute them are based on the template 
fitting of observed spectral energy distributions (SED) and, thus, require multi-colour 
photometry. Here, we demonstrate that, in cases of widely used optical and near- 
infrared filters, if -corrections can be precisely approximated as two-dimensional low- 
order polynomials of only two parameters: redshift and one observed colour. With this 
minimalist approach, we present the polynomial fitting functions for if-corrections 
in SDSS ugriz, UKIRT WFCAM YJHK, Johnson-Cousins UBVR C I C , and 2MASS 
JHK S bands for galaxies at redshifts Z < 0.5 based on empirically-computed values 
obtained by fitting combined optical-NIR SEDs of a set of 10 5 galaxies constructed 
from SDSS DR7 and UKIDSS DR5 photometry using the Virtual Observatory. For 
luminous red galaxies we provide if-corrections as functions of their redshifts only. 
In two filters, g and r, we validate our solutions by computing if -corrections directly 
from SDSS DR7 spectra. We also present a if-corrections calculator, a web-based 
service for computing if-corrections on-line. 

Key words: galaxies: (classification, colours, luminosities, masses, radii, etc.) - galax- 
ies: photometry - galaxies: evolution - galaxies: stellar content - galaxies: fundamental 
parameters 



1 INTRODUCTION 

Extragalactic studies usually require comparison between 
photometric data for different galaxy samples, in particu- 
lar, comparing measurements obtained for distant galaxies 
to the local Universe, where properties of galaxies are stud- 
ied in a much greater detail. Generally, any differences in 
observable parameters arise from: (1) astrophysical prop- 
erties of galaxies and (2) observational biases. The former 
ones include the galaxy evolution effects: due to the light 
travel time, we see distant galaxies as they were looking 
several Gyr ago, so the evolution during the last period of 
their lifetime simply cannot be observed. On the other hand, 
observational biases arise from the process how the observa- 
tions are carried out and, therefore, change drastically from 
one facility to another, even assuming the data are perfectly 
reduced and calibrated. These include effects of aperture, 
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spatial resolution, and a family of effects connected to the 
photometric bandpasses. 

Photometric data often originate from different obser- 
vational studies exploiting different photometric systems 
and, thus require col our transformations to be applied (e.g. 
iFukudta et al.lll995h . But even if the data are obtained in 
the course of one given project, a galaxy sample may contain 
objects at different redshifts. 

In a broad wavelength range, from ultra-violet to 
near-infrared, the spectral energy distributions (SED) of 
non-active galaxies are mostly determined by their stel- 
lar population properties, i.e. age and chemical compo- 
sition, and effects of internal dust attenuation increas- 
ing dramat i cally at short wavelengths ([Calzetti et al 1 ll994l : 
iFitzpatrickl [1999). Stellar population SEDs are very far 
from flat distributions and exhibit prominent features (e.g . 
iFioc Rocca-Volmerangel Il997l : iBruzual fc Charlotl 120031 ) . 
At the same time, redshift ing the galaxy spectrum is equiv- 
alent to shifting the corresponding filter transmission curve. 
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This explains the difference in fluxes in the same band- 
pass for two hypothetical galaxies having exactly identical 
SEDs but being at different re dshifts. Historically, th is dif- 
ference is called X- correct ion (|Oke fc Sandagdll968h . The 
if -correctio n formalism i s presented in detail and thorough ly 
discussed in lHogg et all (l2002h ; IBlanton fc Roweisl (|2QQ7h . 

Nowadays, in the era of large wide-area photo- 
metric and spectroscopic extragalactic surveys, the pre- 
cise, fast and simple computation of iT-corrections 
has become a crucial point for the successful as- 
trophysical interpretation of data. Several approaches 
were presented in t he literature (IFukugita et aD fl 995; 
Mannucci et al.ll200ll; IBlanton Roweisll2007l : iRoche et al.1 
2009h . IBlanton fe RoweisT ([20071 ) provide a software package 
to compute iT-corrections for datasets in any photometric 
system. However, since their method is based on the SED fit- 
ting technique, the results critically d epend on the ayailabil - 
ity o f multi-colour p h otom etric data. IFukugita et al.l ([l995) 
and iMannucci et al.l ([200 lh provide only qualitative depen- 
dence of if- correct ions on redshift and galaxy morphological 
type; the latter is very difficult to assess in an automatic way 
and these methods therefore require availability of original 
galaxy images in addition to photometric measurements. 

The aim of our study is to explore the parameter space 
of typical observed galaxy properties and to provide sim- 
ple and precise analytical approximations of if-corrections 
in widely used optical and NIR photometric bandpasses, 
based on the minimal set of observables. To achieve this goal, 
we exploit a large homogeneous database of optical-to-NIR 
galaxy SEDs compiled from modern wide-area photometric 
surveys. In the next section, we describe our galaxy sample, 
details on the computation of if- correct ions using PEGASE.2 
stellar population models (|Fioc &; Rocca-Volmerange1ll997l ) 
and comparison of obtain ed values with those com puted us- 
ing the kcorrect code (|Blanton Roweid [2007V In Sec- 
tion 3, we describe analytical approximations, and the vali- 
dation of our results using spectral-based if-corrections. In 
Section 4, we compare our results with the literature and 
briefly discuss some astrophysical interpretation of our tech- 
nique. Appendices provide tables with coefficients of best- 
fitting polynomials and present a "K-corrections calculator" 



2 EMPIRICAL COMPUTATION OF 
if-CORRECTIONS 

2.1 Galaxy sample 

We compute if-corrections using a large sample of optical- 
to-NIR SEDs of nearby galaxies constructed using Virtual 
Observatory technologies to retrieve and combine photo- 
metric measurements fro m Sloan Digital Sky Su rvey Data 
Release 7 (SDSS DR7, lAbazaiian et al.1 l2009h and the 
UKI RT Infrared Deep Sur vey Data Release 5 (UKIDSS 
DR5, lLawrence et aljl2007b . Comprehensive description of 
this multi-colour photometric catalogue will be provided in 
a separate paper (Chilingarian et al. in prep), here we give 
a brief summary essential for understanding the empirical 
if-correction computations. 

We constructed a sample of galaxies excluding broad- 
line active galactic nuclei (AGN) by performing the spatial 



cross-matching of the SDSS DR7 spectroscopic sample in 
stripes 9 to 16 with the Large Area Survey catalogue of 
UKIDSS DR5. We selected 190,275 galaxies having spec- 
troscopic redshifts in a range 0.03 < Z < 0.6 provided by 
SDSS DR7 using the SDSS CASJobs ServiccQ. The spatial 
cross-identification with UKIDSS DR5 with a search radius 
of 3 arcsec selecting the best positional matches in case of 
multiple objects within this radii resulted in 170,533 objects, 
87,161 of which were detected in all four UKIDSS Large Area 
Survey photometric bands (Y, J, if, and K). 

In order to compute if-corrections, the photometric 
measurements from the two data sources have to be ho- 
mogeneous. We use SDSS fibre magnitudes corresponding 
to d — 3 arcsec circular apertures (fiber Mags), and com- 
puted the corresponding 3 arcsec aperture magnitudes for 
UKIDSS objects by linearly interpolating between the val- 
ues provided for three apertur es (2.0, 2.8, and 5. 7 arcsec) 
applying zero-point corrections ([Hewett et al.ll2006h for con- 
verting UKIDSS magnitudes from the Vega into the AB sys- 
tem. 

The lower redshift limit, Z — 0.03, is selected in order 
to minimise the aperture effects: at a distance of 120 Mpc, 
corresponding to this redshift, the 3 arcsec aperture encloses 
about 1.75 kpc, i.e. significant part of the bulge even for gi- 
ant galaxies, thus stellar populations in galactic nuclei would 
not dominate the light. Beyond the selected upper redshift 
limit, Z — 0.6, the fraction of normal galaxies in SDSS sig- 
nificantly decreases, because of the magnitude-limited selec- 
tion of SDSS spectroscopic targets, and at the same time 
the quality of absorption-line spectra becomes quite poor. 

We use magnitudes in 3 arcsec apertures and not the 
Petrosian magnitudes to be able to compare them directly 
with the SDSS DR7 spectra obtained within the same aper- 
tures. The median 3 arcsec aperture magnitude uncertain- 
ties are 0.01 mag and better for g, r, z, Y, J, ii, and K, 
0.017 mag for z, and 0.07 mag for the u band respectively. 
All magnitudes are co rrected for the foreg round Galactic 
extinction according to ISchlegel et al.l ([ 19981 ) . 

There are important systematic offsets of unknown ori- 
gin between SDSS DR7 fibre magnitudes and UKIDSS DR5 
photometry in 3 arcsec apertures. We fit 3rd order poly- 
nomial using 5 colours starting from r — i and redder ex- 
cept z — Y and compute the offset between the expected 
z — Y value from the best-fitting polynomial and the ob- 
served one. The offset has a mean value of 0.22 mag and a 
standard deviation of 0.13 mag independent from other pa- 
rameters (e.g. observed colours and a redshift). We therefore 
subtract it from all UKIDSS magnitudes. This effect is illus- 
trated in Fig. [1] In the top panel, we display the combined 
SDSS-UKIDSS SED of some galaxy and its best-matching 
PEGASE.2 template obtained by the iT-correction determi- 
nation procedure described below. The middle and bottom 
panels show the residuals between the observed SED and its 
model, and colours in consequent spectral bands with the 
best-fitting 3rd order polynomial. From the middle panel, it 
is clear that the NIR part of the SED is offset by a constant 
value from the optical part, which is evident as the measured 
z — Y colour strongly deviates from the fitting polynomial. 
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Figure 1. Top panel presents the optical-to-NIR SED of 
SDSS J155023. 03-000023. 8 in a 3-arcsec aperture as provided in 
the catalogue (black diamonds), its NIR part empirically cor- 
rected (blue diamonds), its best-matching PEGASE.2 template 
used for K-correction computation (solid red line). Fitting resid- 
uals are shown in middle panel (same symbols). Bottom panel 
displays colours as a function of wavelength (black crosses), the 
best-fitting 3rd order polynomial for 5 colours (r — i, i — z, Y — J, 
J — H, and H — K) is shown as a solid green line; the corrected 
value of z — Y is denoted by the blue cross. 



2.2 Computation of if-corrections 

We used two approaches to comp ute K- correct ions: ( 1) the 
kcorrect software package by iBlanton & Roweid (|2007h 
and (2) PEGASE.2-based computations described hereafter. 
The latter technique allows us to roughly estimate mean 
stellar population properties as well as internal extinction in 
galaxies. The analytical fitting is performed for both meth- 
ods independently. 

We first ran the kcorrect software package to compute 
the K- correct ions for our sample of galaxies. This package 
is based on a mathematical algorithm, namely non-negative 
matrix factorisation, which creates model-based template 
sets. The initial set of hundreds of templates is reduced to 
a basis of five, which in principle can be used to interpret 
the galaxy SED in terms of stellar populations. The linear 
combination of these templates is then fitted into the set 
of broadband fluxes available for each galaxy to derive the 
i^T- correct ions in all bands. 

Second, we computed a grid of Simple Stellar Popula- 
tions (SSP) using the PEGASE.2 e volutionary synthesis code 
(|Fioc &: Rocca-Volmerangell 19971 ) for a set of 75 ages nearly 
logarithmically spaced between 25 Myr and 16.5 Gyr and 10 
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Figure 2. Comparison of empirically calculated X-corrections 
obtained by fitting the photometric data using PEGASE.2 SSP 
models (kp^) and by the KCORRECT (kBR07) code. Each panel 
displays the difference between the two approaches as a function 
of redshift. Solid lines denote the median differences and their 
standard deviations are shown with dotted lines. 

metallicities between —2.5 < [Fe/H] < +1.0 dex. Such a grid 
was computed separately for reds hifts between < Z < 0.6 
with a step of 0.05. We apply the iFitzpatrickl (|1999h extinc- 
tion law for each of 750 SSPs at each redshift varying the 
Ay between and 2.25 mag with a step of 0.15 mag ending 
up with 11250 template SEDs per redshift. 

In order to compute if- correct ions for a given galaxy, 
we perform a linear interpolation of the SSP grid to its red- 
shift, then pick up the best-matching template SED in terms 
of x 2 normalising both data and templates by the mean 
fluxes in all filters. Since our photometric uncertainties are 
quite small, this approach would not result in significant 
biases. Once the best-matching template has been found, 
the K- correct ions in all photometric bands are computed as 
K f (Z) = -2.51og(F(0)/F(Z)), where F(Z) and F(0) are 
fluxes in a given filter at redshift Z and in the restframe. 
We provide an example of an observed SED and its best- 
matching template in the upper panel of Fig [1] 

Since we can calculate model fluxes in any photomet- 
ric band, we also used the same technique to compute 
if-corrections in Johnson- Cousins UBVR C I C and 2MASS 
JHK S bands. 

We compare the values computed in this fashion with 
the value derived by the kcorrect code from the same pho- 
tometric dataset. The comparison for all 9 SDSS-UKIDSS 
bands is provided in Fig [2) 

In general, the results obtained by using the two ap- 
proaches in ugrizY JHK bands are quite similar. However, 
in certain spectral bands, some statistically significant dif- 
ferences are evident. The worst situation is observed in the 
SDSS u band mainly for two reasons: relatively poor qual- 
ity of the u band photometry especially for objects at higher 
redshifts and very high sensitivity of UV colours to even low 
mass fractions of recently formed stars, which affect fluxes at 
longer wavelengths much weaker. Therefore, the SSP fitting 
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approach presumably should not work well. At the same 
time, the linear combination of 5 templates used by the 
KCORRECT package may also produce significant biases due 
to age and/or metallicity mismatch between the templates 
and real galaxies. Finally, we provide the K- correct ions for 
the u band, emphasising that since no independent verifi- 
cation is possible in our case, one has to use these results 
with caution. The same statement applies to K band K- 
corrections, although the results of the two approaches well 
match each other, because the computation relies on the ex- 
trapolation of galaxy SED in the NIR part, where stellar 
population models and, correspondingly, the template spec- 
tra are of much lower quality than in the optical wavelength 
domain. 

Computations in the r, i, and J bands agree remark- 
ably well except the high-redshift (Z > 0.5) end of the J 
band. There are some systematic discrepancies between the 
two techniques in the g, z, Y, and H bands. They may orig- 
inate from the fact that PEGASE.2 SSPs are built using the 
theoretical stellar library, shown to introduce colour differ- 
ences between synthetic spectra and observed ones at least 
in the SDSS photometric system ([Maraston et al.l l2009). We 
will analyse the g and r band results below using direct 
spectral-based iT-corrections, while for the remaining band- 
passes no independent test can be performed since no large 
samples of galaxy spectra are available in those wavelength 
domains. However, we note that the discrepancies are an or- 
der of 0.05 mag, hence, both approaches may be used in 
the photometric studies. Therefore, we will proceed with 
the rest of our analysis using both techniques of the in- 
correct ion computation, addressing them as BR07 and SSP 
for the KCORRECT and pegase.2 SSP-based approaches re- 
spectively. 



3 RESULTS 

3.1 Analytical approximations 

We observe a large scatter of K- correct ions as functions of 
redshift reaching 2 mag in all SDSS-UKIDSS bands except 
H and K. However, exploring the data with the topcat soft- 
war^E we found that adding just one observed colour as a 
second parameter and approximating the if-correction as a 
surface in the three-dimensional space, significantly reduces 
the residual scatter bringing it to the order of iT-correction 
computation uncertainties. We fit if -correction values in ev- 
ery filter q as a polynomial surface of a form: 

N z N c 

K q (Z,m fl -m h ) =^2^2a x , y Z x (m fl -m h ) y , (1) 

cc = y = 

where ^x,y are polynomial coefficients, Z is a spectroscopic 
redshift, mf 1 and m/ 2 are observed magnitudes in filters 
fi and fi chosen for every filter q, Nz and N c are empiri- 
cally selected polynomial powers in the redshift and colour 
dimensions respectively. Given that if-corrections are zero 
by definition at Z = 0, no constant term is needed and all 

CL0,y — 0. 
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Since both if-correction computation techniques used 
in our study are based on the stellar population models, and 
our analytical approximations exploit simple polynomial fit- 
ting technique without any clipping of outliers, we have to 
exclude strong active galactic nuclei (AGN) and quasars, 
as well as SDSS targeting artifacts (e.g. aircrafts, satellites, 
minor planets) and objects with wrong redshift determina- 
tions, which may affect our best-fitting solutions. We have 
fitted all 190,275 SDSS D R7 spectra using the nbu rsts full 
spectral fitting technique (| Chilingarian et aD2007bl lah in the 
rest frame wavelength range between 3900 and 6750 A, and 
selected for our further analysis only 164,108 objects having 
reduced x 2 /DOF < 0.9 (median X 2 /DOF = 0.67, which is 
smaller than unity because of slight spectral oversampling of 
SDSS data). The interpretation of the spectral fitting results 
will be provided together with the presentation of a spec- 
trophotometric catalogue in Chilingarian et al. (in prep). 

We computed and fitted if-corrections in this fashion 
for a sample of 164,108 SDSS DR7 galaxies with well-fitted 
spectra using only SDSS 5-bands photometry, and to the 
merged SDSS-UKIDSS sample containing photometric in- 
formation in all 9 bands for 74,254 galaxies. While the re- 
sults from the two approaches remain nearly statistically 
identical in u, g, and r bands, the difference becomes signif- 
icant in i and especially in z. Two distinct sequences become 
clear in the redshift vs if-correction plots in i and z in case 
of the 5-bands based computation. Only one of the two se- 
quences remains in each case if full 9-band SEDs are used, 
suggesting that the second sequence is created by wrong stel- 
lar population templates picked up from the template grid. 
This clearly demonstrates the importance of the NIR part to 
compute K- correct ions if using a fully empirical approach. 
Consequently, r.m.s of the polynomial surface fitting resid- 
uals in case of 9-bands based K- corrections in the z band is 
four times smaller compared to the 5 bands. 

In Fig [3HH we display the analytical approximations 
of if-corrections in 9 bands for the values computed us- 
ing PEGASE.2-based matching and the KCORRECT code. Up- 
per panels demonstrate computed values as a function of 
redshift and an observed colour, colour-coded according to 
the scale bars presented in the plots. The lower panels dis- 
play the mean residuals and their r.m.s. as a function of 
redshift. Since we use a fixed grid of pegase.2 templates 
and do not perform interpolation on the age and metal- 
licity axes, the individual template age sequences become 
visible at high redshifts range in some filters (e.g. Y). How- 
ever, the computational errors due to this discretisation do 
not exceed 0.03 mag. Black lines in the plots denote the 
behaviour of if -corrections for galaxies having fixed rest- 
frame (i.e. if-corrected) g — r colours. Solid lines are for 
0.73 < g-r < 0.81 mag, dashed for 0.58 < g-r < 0.70 mag, 
dot-dashed for 0.4 < g — r < 0.6 mag, and triple-dot-dashed 
for g — r < 0.15 mag, roughly corresponding to luminous red 
galaxies (LRG), early- type spirals (S-early), late- type spi- 
rals (S-late), and actively star-forming galaxies (SF). These 
lines are constructed by connecting median values within 
0.03- wide redshift bins in the colour bins given above. We 
notice that the i^-band if-corrections computed with the 
KCORRECT code at higher redshifts turn to be higher than 
expected for LRGs defined by g — r rest frame colours, which 
is probably indicative of a template mismatch. 
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Figure 3. X-corrections in 9 photometric SDSS-UKIDSS bands computed using PEGASE.2 SSP matching and the residuals of the 
analytical fitting. Upper panel for each filter presents computed X-corrections vs redshift with a colour-coded observed colour used to 
perform the fitting. Sequences of galaxies with constant restframe g — r colours roughly corresponding to different morphological types 
are overplotted (see the text for details). Lower panels display the residual between the analytical approximation and measured values 
(bold black line) and r.m.s. of the fitting residuals (dotted lines). 



The 4 morphological types are well separated in the u, 
g, and r bands; only very blue objects exhibit a distinguished 
behaviour in z, while in NIR bands the i^- correct ions become 
nearly independent of a restframe colour of a galaxy (or its 
morphological type). 

Some tables containing the coefficients of polynomial 
approximations are provided in Appendix A. Additional 
tables for if- correct ion approximations using different ob- 
served colours are provided on the "K-corrections calcula- 
tor" web-site described in Appendix B. 

Classical Johnson-Cousins photometric system is still 
widely used for extrag alactic research, as w ell as NIR bands 
of the 2MASS survey (jSkrutskie et al.ll2QQ6h . therefore com- 
putation and analytical approximation of K- correct ions in 



these bands are of a great practical importance. However, we 
do not have photometric measurements for galaxies in our 
sample obtained in these bands. We used th e photometric 
transformations denned in Ijordi et al.l (|20 06) for stars and 
those provided on the web-pages of SDSS in order to con- 
vert available SDSS ugriz magnitudes into Johnson-Cousins 
UBVRcIc, and UKID SS YJHK to 2MAS S JHK S trans- 
formations presented in lrlewett et al.l (|2QQ6h . Then, we used 
these magnitudes to fit if- correct ions in the UBVR C I C and 
2MASS JHK S bands as functions of "computed" observed 
colours denned in the same photometric system. The coef- 
ficients of these polynomial approximations are available at 
the "K-corrections calculator" web-site. We note that the 
approximations for the Johnson-Cousins and 2MASS niters 
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Figure 4. The same as Fig [3] but the X-corrections are computed using the KCORRECT code. 



are provided as functions of colours in the Vega system, 
whereas for SDSS and UKIRT WFCAM bands the colours 
are expressed in AB magnitudes. 

3.2 Validation using SDSS DR7 spectra 

We fetched optical SDSS spectra in the wavelength range 
between 3800 and 9200 A for all galaxies from our sample, 
therefore we were able to perform independent direct ver- 
ification of the analytical approximations of if- correct ions 
presented above. The available spectral coverage allows us 
to compute fluxes directly from the spectra using both real 
and redshifted filter transmission curves, i.e. to obtain real 
rest frame fluxes, for g, r, and i filters. The highest redshifts, 
where reliable computation is still possible for the r and i 
bands are ~ 0.28 and ~ 0.08 respectively due to the upper 
wavelength limit of the SDSS spectral c overage. We use a 
similar technique to iRoche et al.l (|2009 ) for the computa- 



tion of spectral-based i^- correct ions, with two main differ- 
ences: (1) we redshift the filter transmission curves instead 
of blueshifting the spectra in order to prevent interpolation 
of fluxes and smoothing defects in the spectra; (2) in case 
of missing data within the wavelength range, we interpolate 
fluxes linearly or extrapolate them using constant level in 
F\ if the red tail of the redshifted filter transmission curve 
goes beyond 9200 A. The maximal adopted truncation of the 
filter transmission curve was at a level of 1 per cent of the 
maximal transmission. 

In Fig we present the differences between spectral- 
based if-corrections in g and r bands and the analytical 
approximations of those computed empirically from SDSS 
fiberMag using pegase.2 SSP templates and the KCORRECT 
code as functions of redshift and observed g — r colour. In the 
r band, the agreement between spectral-based if -correct ions 
and analytical fitting derived above is as good as 0.02 mag. 
For the low-redshift (Z < 0.3) part of the sample, the situ- 
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at ion is similar in the g band, although the systematic dif- 
ference reaches 0.06 mag. However, for the higher redshifts 
(0.3 < Z < 0.5), the discrepancies become as large as 0.15- 
0.2 mag. This inconsistency possibly originates from under- 
estimated g synthetic magnitudes due to a very low signal 
in the blue part of the spectra for objects having Z > 0.3. 
Formally computed photometric errors of synthetic g mag- 
nitudes in this redshift range reach 0.2 mag. 



4 DISCUSSION AND CONCLUSIONS 
4.1 Comparison with literature 

We compare if- correct ions com puted in this study wit h the 
results described in literature. IFukugita et al.l (|l995h pre- 
sented if- correct ions of the SDSS photometric system in 
their fig. 20. Since no data are provided in the numeri- 
cal form, we are able to compare the results only quali- 
tatively. The behaviour of the values as functions of red- 
shifts agree well at Z < 0.5 , differ ent morphological types of 
galaxies in IFukugita et~aT1 (|l995l ) are related to the recon- 
st ructed rest fram e colo urs. The elliptical galaxy template 
of IFukugita et al.1 (|l995l ) behaves similarly to the sequences 
denoted as "LRG" in Fig. [3HI] corresponding to luminous 
red gal axies. We notice that the filter transmission curves 
used in IFukugita etHI (jl995h are somewhat different from 
those obtained from the telescope and presented at the SDSS 
DR7 web pages, especially for the g and r bands, explaining 
why the exact match between the two approaches cannot be 
achi eved. 

iMannucci et al.1 (|200lh present in their fig. 7 NIR K- 
corrections computed from the templates spectra of galaxies 
having different morphological types computed with the pe- 
GASE code. Although their JHK bands are slightly different 
from those of the UKIDSS, the behaviour of if- correct ions as 
functions of redshift is very similar. Authors mention that 
if-corrections in H and K are virtually insensitive to the 
galaxy morphological types, which we observe as their very 
weak dependence on galaxy colours. 

We compare our re sults with the spe ctral-based K- 
corrections presented in iRoche et al.l (|2009h . The authors 
deal only with red sequence galaxies, therefore we are able 
to test only a small although very important part of the pa- 
rameter sp ace. The values of if -corrections provided in ta- 
bles 1-2 of lRoche etHI (l2009h are overplotted in Fig [5] The 
agreement is remarkably good, which is, however, expected 
because our spectral-based if-corrections computed in the 
same manner agree well with the best-fitting solutions. 



4.2 Recovered restframe magnitudes and LRGs 

It is well known that the restframe colours correlate with the 
galaxy morphology, and there is a pronoun ced bimodality 
in the colour distribution of galaxies (see e.g . IStrateva et al.l 
l200ll : iBaldrv et al.ll2004l : iBalogh et al.ll2004h . The so-called 
"red sequence" includes elliptical and lenticular galaxies 
with some fraction of early-type spirals demonstrating a 
lack of ongoing star formation, whereas the "blue cloud" 
contains actively star- forming objects represented mostly 
by late- type spiral and irregular galaxies. Some morpholog- 
ically classified early-type galaxies indeed sit below the red 
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Figure 5. Differences between analytically approximated and di- 
rectly measured spectral-based X-corrections in the g (top) and 
r (bottom) bands. Upper panels display computed values of K- 
corrections. Dashed lines denote spectral-based if -co rrections for 
early- type galaxies presented in I Roche et al.l (|2Q09h . Two bot- 
tom panels for each spectral band display the differences between 
spectral-based K-corrections and those computed using PEGASE.2 
SSP matching and the KCORRECT code, respectively. 
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Table 1. Coefficients of the best-fitting polynomials for the 
redshift-dependence of if -correct ions for luminous red galaxies. 
Initial if -correction values were computed from the PEGASE.2 
SSPs. 





Z 1 


Z 2 


Z 3 


Z 4 


Z 5 


K u 


5.93938 


-30.5247 


179.473 


-380.488 


282.011 


Kg 


2.61617 


-4.44391 


93.0132 


-284.582 


252.245 


K r 


0.312233 


14.3325 


-68.2493 


136.254 


-87.3360 


Ki 


0.234538 


14.3162 


-97.2754 


246.775 


-207.028 


K z 


0.897075 


3.60112 


-34.7890 


93.0266 


-79.5246 


K Y 


0.402992 


5.30858 


-18.3172 


17.6760 


-0.31400 


Kj 


-0.076704 


-4.48411 


27.1585 


-45.6481 


22.6928 


K H 


0.382926 


-1.81590 


-13.1657 


57.5486 


-59.0677 


K k 


-1.75997 


5.48023 


-56.4175 


175.939 


-160.754 



Figure 6. Colour-magnitude diagram presenting data for 74,254 
galaxies using the SDSS DR7 and UKIDSS DR5 Petrosian mag- 
nitudes in the g, r, and H bands respectively If -corrected us- 
ing the analytical approximations presented in this study. Black 
and blue contours are for the PEGASE.2-based and KCORRECT de- 
rived values respectively. Levels of the two-dimensional density 
plot correspond to the powers of two from 2 (outermost) to 1024 
(innermost) galaxies per bin of 0.25 magxO.025 mag. 



sequence in the "green valley" or even "blue cloud" regions. 
In Fig [6] we provide the colour-magnitude plot of 74,254 
galaxies from SDSS DR7 and UKIDSS, where their total 
Petrosian magnitudes were if-corrected using the presented 
analytical approximations. 

Some of the objects residi ng below the red se- 
quence referred as E+A galaxies (IPressler Gurml Il983l ; 
ICouch fe Sharpleslll987l : IPressler et alJ Il999h often exhibit 
weak if any emission lines in their spectra ruling out major 
star formation. However, Balmer absorptions are remark- 
ably deep, which is a good evidence for the presence of 
young stars that was confirmed by the recent study wher e 
authors used full-spectral fitting (IChilingarian et alJl2009bh . 
Despite the early- type morphologies of E+A, these galax- 
ies are indistinguishable from spirals if one uses only in- 
tegrated colours and therefore, colour transformations and 
if- correct ions behave very similarly for these two morpho- 
logically different classes of galaxies. Hence, the rest frame 
colours are not tightly connected to the galaxy morphol- 
ogy, but to the internal properties such as ongoing star for- 
mation and, therefore, they are specific of stellar popula- 
tions. Another example would be dwarf elliiptical (dE) and 
compact elliptical (cE) galaxies. The former ones usually 
have intermediate- age stellar popu lations with metallicities 
lower than those o f LRGs (see e.g. IChilingarian et aL l l2008l : 
IChilingarianl 12009) , therefore their global optical colours are 
bluer, and their if-corrections in the corresponding bands 
are lower than those of more massive early- type galaxies, 
which results in quite a strong tilt of the red sequence at 
Mh,ab > —20 mag (se e Fig [6]) reproducin g its behaviour 
for nearby dE galaxies ([janz &: Lisked l2009h . On the other 
hand, cE galaxies with luminosities similar to dEs, originate 
from the tidal stripping of more massive early-type progeni- 
tors and, therefore, exhibit old metal-rich stellar populations 
(Chi lingarian et aUl2009al ). For this reason, cEs have signif- 
icantly redder colours compared to dEs and, therefore, have 
behaviour of if- correct ions very similar to LRGs. 

There is a weak dependence of the colour of red se- 



Table 2. Same as in Table [T] but the KCORRECT computed K- 
corrections. 

Z 1 Z 2 Z 3 Z 4 Z 5 



K u 


4.20000 


-24.5015 


229.149 


-574.272 


434.901 


Kg 


2.17470 


10.3810 


1.49141 


-76.6656 


88.6641 


K r 


0.710579 


10.1949 


-57.0378 


133.141 


-99.9271 


Ki 


0.702681 


4.27115 


-37.2060 


112.054 


-105.976 


K z 


0.643953 


-1.88400 


13.6952 


-35.0960 


29.9249 


K Y 


-0.245996 


21.8772 


-137.019 


322.051 


-257.136 


Kj 


0.106358 


-5.06024 


18.4707 


-3.73196 


-23.9595 


K H 


0.268479 


3.03488 


-35.8994 


98.6524 


-83.9401 


K K 


-2.80894 


15.6923 


-96.8401 


256.235 


-220.691 



quence galaxies on the galaxy luminosity reflecting, in par- 
ticular, a mass-metallicity relation of early-type galaxies. 
Massive (and luminous) red sequence galaxies compared to 
lower-mass systems contain more metal-rich stellar popula- 
tions having intrinsically redder colours than the metal-poor 
ones. 

The special case of LRGs is very important for under- 
standing galaxy formation and evolution, therefore we pro- 
vide specific approximations for these objects. Since the in- 
trinsic spread of LRGs restframe colours is very low, about 
0.04 mag, their if -corrections can be well approximated as 
a function of a single parameter, the redshift. We selected a 
sample of LRGs based on their restframe colours, and fitted 
the if- correct ions as polynomial functions of their redshifts 
in all 9 photometric bands. The coefficients of best-fitting 
5th order polynomials with no constant term are presented 
in Table [T]and Table [2] for the pegase.2 SSP-based and 
KCORRECT-computed values correspondingly. Every row cor- 
responds to a given photometric band and contains coeffi- 
cients from the 1st to the 5th power of redshift. 

Rest-frame colours of LRGs provide a natural test for 
the quality of if- correct ions. The vast majority of stars 
in these objects is believed to form on short timescales in 
the early Universe and then evolve passively. The redshift 
z = 0.5 corresponds to the lookback time about 5 Gyr. Thus, 
if one assumes LRGs to be as old as 12 Gyr in the local Uni- 
verse and contain no younger stars, then at z — 0.5 their 
restframe optical colours would be slightly bluer with the 
differences A(u — r) w 0.25 mag, A(g — r) w 0.08 mag, 
and redder colours different by less than 0.03 mag (values 
estimated from the colour evolution of the solar metallicity 
PEGASE.2 SSP). In Fig [7] we present the behaviour of recov- 
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Figure 7. Dependence of recovered restframe colours of LRGs on 
redshift. Solid and dashed lines are for the PEGASE.2 SSP-based 
and KCORRECT computed if -corrections respectively. 

ered restframe colours of LRGs as functions of redshift. All 
presented colours have SDSS r as one of the bands, because 
the best consistency of both if -correction computation al- 
gorithms is reached in this band. 

The bluest u — r colour exhibits significant changes over 
redshift exceeding by a factor of 3 the expectations from the 
passive evolution of a SSP. This, however, can be explained 
by a "tail" of the star formation history since even small 
mass fractions of intermediate mass stars strongly affect the 
u band photometry. On the other hand, we cannot exclude 
the template mismatch between the model and real galaxies 
to be partly responsible for this effect. 

Surprisingly, the g — r colour does not evolve at all if 
we use the KCORRECT values of g-band if -correct ions, but 
evolves slightly stronger than expected from the SSP evolu- 
tion when using the PEGASE.2 based values, being somewhat 
consistent with the u — r colour behaviour. 

While r — i and r — z colours behave similarly in both 
approaches exhibiting, as expected, virtually no evolution, 
in NIR bands the situation is different. Both techniques are 
consistent in the J band, although demonstrating rather un- 
expected evolution by about 0.08 mag, whereas in H and if 
they are significantly different (as also can be seen from 
the corresponding panels in Fig. [2}. The pegase.2 SSP- 



based if-band if-corrections at Z > 0.3 look slightly un- 
derestimated, conversely to the KCORRECT-computed ones, 
which are overestimated. In the H band, both methods seem 
to underestimate if-corrections, but SSP-based approach is 
stronger affected. 

The decisive answer to the questions about the K- 
correction computation in NIR bands will be given only 
when the next generation NIR stellar population models 
based on empirical stellar libraries become available. The 
H and if band solutions presented in this paper have to 
be used with caution keeping in mind that some systematic 
errors may be introduced to the final results. 

4.3 Summary 

We present polynomial approximations of if- correct ions in 
commonly used optical (SDSS ugriz and Johnson-Cousins 
UBVRJc) and near-infrared (UKIRT WFCAM YJHK 
and 2MASS JHK S ) broad-band filters as functions of red- 
shift and observed colours for galaxies at redshifts Z < 0.5. 
The traditional if-correction computation techniques based 
on the SED fitting require multi-colour photometry, which 
is not always available. Our approach allows one to com- 
pute restframe galaxy magnitudes of the same quality using 
a minimal set of observables including only two photometric 
points and a redshift. For luminous red galaxies, we pro- 
vide the fitting solutions of if-corrections as 1-dimensional 
polynomial functions of a redshift. 
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Table Al. The coefficients a x , y of the two-dimensional polyno- 
mial approximation of if -corrections K u (Z,u — r) in the SDSS 
u band (see Equation 1). In order to derive a value for K u , one 
needs to sum the polynomial terms of a form Z x and (u — r) y 
from the table heading and its first column multiplied by the 
coefficients in corresponding table cells. 





(u - r)° 


(u — r) 1 


(u — r) 2 


(u — r) 3 
















z 1 


1.63349 


2.24658 


0.141845 


-0.13441 


z 2 


-71.84 


20.4939 


-3.82771 


0.789867 


z 3 


257.509 


-42.3042 


-4.05721 





z 4 


-308.573 


63.0036 








z 5 


42.8572 











Table A2. Same 


as in Table lAll but for the SDSS g 


function of redshift and the g 


— r colour. 






(g-r)° 


(g-r) 1 


{g-r? 


(g-r) 3 


Z° 














z 1 


-0.900332 


3.97338 


0.774394 


-1.09389 


z 2 


3.65877 


-8.04213 


11.0321 


0.781176 


z 3 


-16.7457 


-31.1241 


-17.5553 





z 4 


87.3565 


71.5801 








z 5 


-123.671 












in the Johnson-Cousins UBVR C I C bands. The coefficients 
for the 2MASS JHK S bands as well as for the other colour 
combinations of SDSS and UKIDSS bands are provided on- 
line through the "K-corrections calculator" service. 



APPENDIX A: 2-DIMENSIONAL 
APPROXIMATIONS OF if-CORRECTIONS 

Here we provide tables containing coefficients of the ana- 
lytical approximations of if-corrections. Columns in every 
table contain colour terms of the fitting solution, from con- 
stant to the 3rd order respectively. Rows are for the red- 
shift dependence. The maximal degree of the polynomial 
surface is set to 5 or 7 (for the H band), therefore the 
coefficients for the higher degrees, in both colour and red- 
shift, are set to zero. Tables [ATHA91 and IA1QHA18I present 
the coefficients for approximations of PEGASE.2-based and 
KCORRECT-computed if -corrections in the SDSS ugriz and 
UKIDSS YJHK bands. Tables [A19HA231 contain the coeffi- 
cients for approximations of PEGASE.2-based if -corrections 



Table A3. Same as in Table lAll but for the SDSS 
function of redshift and the g — r colour. 



band as a 





(g-r)° 


(9-r) 1 


(g-r) 2 


(g-r) 3 
















z 1 


-1.61294 


3.81378 


-3.56114 


2.47133 


z 2 


9.13285 


9.85141 


-5.1432 


-7.02213 


z 3 


-81.8341 


-30.3631 


38.5052 





z 4 


250.732 


-25.0159 








z 5 


-215.377 
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Table A4. Same as in Table [All but for the SDSS 
function of redshift and the g — i colour. 



band as a Table A9. Same as in Table lAll but for the UKIDSS K band as 

a function of redshift and the J — K colour. 





(g-i)° 


(g-i) 1 


(g-i) 2 


(g-i) 3 




(J-Kf 


(J-K) 1 


(J-K) 2 


(J-K) 3 


z° 














Z° 














z 1 


-2.41799 


4.68318 


-3.70678 


1.5155 


z 1 


-2.80374 


4.14968 


1.15579 


-1.94003 


z 2 


11.2598 


5.14198 


-2.64767 


-3.63215 


z 2 


13.4077 


-39.5749 


11.7 


4.7809 


z 3 


-94.7387 


-14.154 


27.2864 





z 3 


-69.7725 


94.0769 


-35.1023 





z 4 


285.775 


-51.6662 








z 4 


157.649 


-44.0291 








z 5 


-222.641 











z 5 


-132.317 












Table A5. Same as in Table [All but for the SDSS z band as a 
function of redshift and the r — z colour. 





(r - z)° 


(r-z) 1 


(r - z) 2 


(r-z) 3 


z° 














z 1 


-1.7252 


3.35566 


0.469411 


0.350873 


z 2 


14.9772 


-23.1956 


-3.32427 


1.78842 


z 3 


-41.1269 


75.2648 


-10.2986 





z 4 


11.3667 


-48.3299 








z 5 


23.5438 












Table A 10. Same as in Table lAll but using the kcorrect code. 





(u - r)° 


(u-r) 1 


(u - r) 2 


(u - r) 3 


z° 














z 1 


8.81624 


-12.0027 


5.57928 


-0.825005 


z 2 


33.0392 


39.7152 


-18.7077 


4.0901 


z 3 


-1223.73 


236.944 


-52.6404 


-3.43017 


z 4 


6304.63 


-659.764 


162.85 


-1.62112 


z 5 


-15288.5 


395.301 


-83.6382 





z 6 


18533.6 


-67.5999 








z 7 


-8719.48 












Table A6. Same as in Table [AT] but for the UKIDSS Y band as 
a function of redshift and the Y — H colour. 





(Y - H)° 


(y-H) 1 


(Y - H) 2 


(Y - if ) 3 


z° 














z 1 


-2.01575 


2.70429 


6.01384 


-5.17119 


z 2 


14.3112 


-13.2354 


-4.03961 


13.2633 


z 3 


-46.4835 


12.7464 


-26.3263 





z 4 


80.2341 


27.9842 








z 5 


-66.1958 












Table All. Same as in Table lALOl but for the SDSS g band as a 
function of redshift and the g — r colour. 

(g-r)° (g-r) 1 (g - r) 2 (g - r) 3 



Z° 
Z 1 

z 2 
z 3 
z 4 
z 5 





-0.962084 
15.6602 

-82.9388 
273.308 

-312.677 





2.2796 
-14.8073 
-49.2478 
131.339 






4.16029 
19.261 
-40.9139 









-3.27579 
4.28022 









Table A7. Same as in Table |AT] but for the UKIDSS J band as 
a function of redshift and the J — K colour. 





(J-Kf 


(J-KY 


(J-K) 2 


(J-K) 3 


z° 














z 1 


-0.765217 


2.43055 


-0.427304 


0.277662 


z 2 


1.59864 


-14.646 


12.0911 


-1.2131 


z 3 


-4.02136 


18.077 


-26.1137 





z 4 


18.5608 


25.2691 








z 5 


-40.3567 












function of redshift and the g 


— r colour. 






(g-rf 


(g-r) 1 


(g-r) 2 


(g - r) 3 


Z° 














z 1 


-0.351251 


2.61848 


-2.99032 


1.59058 


z 2 


1.93312 


16.0682 


-2.16736 


-4.24709 


z 3 


-69.9339 


-49.337 


22.9267 





z 4 


253.373 


12.0421 








z 5 


-235.32 












Table A8. Same as in Table|AT]but for the UKIDSS H band as 
a function of redshift and the H — K colour. 





(H -K)° 


(H-K) 1 


(H-K) 2 


(H-K) 3 


z° 














z 1 


-0.642942 


-1.05192 


-15.5123 


-18.1957 


z 2 


26.3667 


80.0291 


192.688 


179.956 


z 3 


-274.11 


-564.952 


-848.543 


-646.653 


z 4 


1081 


1569.47 


1741.31 


761.264 


z 5 


-1938.48 


-1893.45 


-1488.4 





z 6 


1448.38 


869.396 








z 7 


-249.952 












Table A13. Same as in Table lAlOl but for the SDSS % band as a 
function of redshift and the g — i colour. 





(9-i)° 


(g-i) 1 


(g-^) 2 


(g-^f 


z° 














z 1 


-5.58619 


8.51897 


-3.44266 


0.823935 


z 2 


22.398 


-31.8699 


7.13812 


-2.26748 


z 3 


-16.5911 


40.9251 


6.77963 





z 4 


-12.0117 


-44.6466 








z 5 


21.0947 












12 /. Chilingarian et al 



Table A14. Same as in Table lAlOl but for the SDSS z band as a 
function of redshift and the r — z colour. 





(r - z)° 


(r-z) 1 


(r-z) 2 


(r-z) 3 


z° 














z 1 


-1.426 


3.08833 


-0.726039 


1.06364 


z 2 


2.9386 


-8.48028 


-8.18852 


-1.35281 


z 3 


8.08986 


53.5534 


13.6829 





z 4 


-93.2991 


-77.1975 








z 5 


133.298 












Table A15. Same as in Table [AToI but for the UKIDSS Y band 
as a function of redshift and the Y — H colour. 





(Y - H)° 


(Y-H) 1 


(Y - H) 2 


(Y - if) 3 


z° 














z 1 


-2.62137 


4.7578 


2.28856 


-4.02782 


z 2 


29.4209 


-26.5297 


10.1083 


10.5582 


z 3 


-141.372 


34.7785 


-41.2725 





z 4 


311.42 


26.8742 








z 5 


-264.997 












APPENDIX B: K-CORRECTIONS 
CALCULATOR 

To facilitate the usage of described analytical approx- 
imations and calculation of necessary i^- correct ions for 
user's data, we provide a dedicated web-site entitled "In- 
correct ions calculator"0 offering several ways of operation. 
Firstly, one can determine if- correct ions for objects of inter- 
est one by one by manually filling interactive web form with 
the data available. This is the most straightforward method 
to determine if -corrections and it works in any modern web 
browser. As a second option, the ready-to-use code snip- 
pets written in popular data languages (C, IDL, Python) 
together with the corresponding analytical expressions are 
provided to simplify the integration of if- correct ion func- 
tionality into the user's code. Both ways allow a user to cal- 
culate if-corrections for SDSS ugriz and UKIRT WFCAM 
YJHK (in AB magnitudes), but also for Johnson-Cousins 
UBVRJc and 2MASS JHK S bands (in Vega magnitudes), 
choosing the most convenient colour as an input argument 
from several options. Based on user demands and feedback, 
we will consider the development of a web service for bulk 
calculation of if-corrections and publication of additional 
filter-colour combinations in the calculator. 



3 http://kcor.sai.msu.ru/ 

Table A16. Same as in Table lAlO] but for the UKIDSS J band 
as a function of redshift and the J — K colour. 





(J-K)° 


(J-K) 1 


(J-K) 2 


(J-K) 3 


z° 














z 1 


-0.472236 


2.1536 


0.811858 


-1.87211 


z 2 


-0.107502 


-8.00546 


16.6955 


4.85177 


z 3 


-18.2002 


-27.5709 


-46.9334 





z 4 


111.89 


99.1294 








z 5 


-162.057 












Table A17. Same as in Table [Ami but for the UKIDSS H band 
as a function of redshift and the H — K colour. 





(H — K)° 


(H-K) 1 


(H-K) 2 


(H-K) 3 


z° 














z 1 


0.132484 


1.93083 


-4.7581 


-6.67018 


z 2 


8.61784 


-17.3496 


28.7714 


8.05742 


z 3 


-83.7003 


267.725 


127.193 


109.678 


z 4 


134.398 


-1657.94 


-851.199 


-217.543 


z 5 


567.048 


4005.55 


1071.17 





z 6 


-2000.79 


-3298.45 








z 7 


1697.02 











Table A18. Same 


as in Table lAlOl but for the UKIDSS K 


as a 


function of redshift and the J — K colour. 




(J-K)° 


(J-K) 1 


(J-K) 2 


(J-K) 3 


Z° 














z 1 


-3.1771 


1.66876 


1.45967 


-3.40684 


z 2 


17.9897 


-7.83528 


13.3436 


9.32974 


z 3 


-114.067 


-17.793 


-42.0747 





z 4 


318.424 


70.0829 








z 5 


-299.557 












Table A 19. Same as in Table lAll but for the Johnson U band as 
a function of redshift and the U — R c colour. 





(U - R c )° 


(U-Rc) 1 


(U-R c ) 2 


(U - R c ) 3 


z° 














z 1 


2.84791 


2.31564 


-0.411492 


-0.0362256 


z 2 


-18.8238 


13.2852 


6.74212 


-2.16222 


z 3 


-307.885 


-124.303 


-9.92117 


12.7453 


z 4 


3040.57 


428.811 


-124.492 


-14.3232 


z 5 


-10677.7 


-39.2842 


197.445 





z 6 


16022.4 


-641.309 








z 7 


-8586.18 











Table A20. Same 


as in Table [All but for the Johnson B ban 


a function of redshift and the B 


— R c colour. 






(B - R c )° 


(B-Rc) 1 


(B-R c ) 2 


(B-Rc) 3 


Z° 














z 1 


-1.99412 


3.45377 


0.818214 


-0.630543 


z 2 


15.9592 


-3.99873 


6.44175 


0.828667 


z 3 


-101.876 


-44.4243 


-12.6224 





z 4 


299.29 


86.789 








z 5 


-304.526 












Table A21. Same as in Table [All but for the Johnson V band as 
a function of redshift and the V — I c colour. 





(V-Ic)° 


(V-Ic) 1 


(V-Ic) 2 


(V-Ic) 3 


z° 














z 1 


-1.37734 


-1.3982 


4.76093 


-1.59598 


z 2 


19.0533 


-17.9194 


8.32856 


0.622176 


z 3 


-86.9899 


-13.6809 


-9.25747 





z 4 


305.09 


39.4246 








z 5 


-324.357 












Analytical approximations of K -corrections 



Table A22. Same as in Table [All but for the Cousins R c band 
as a function of redshift and the B — R c colour. 





(B - R c )° 


(B-Rc) 1 


(B - R c ) 2 


(B - R c ) 3 
















z 1 


-2.83216 


4.64989 


-2.86494 


0.90422 


z 2 


4.97464 


5.34587 


0.408024 


-2.47204 


z 3 


-57.3361 


-30.3302 


18.4741 





z 4 


224.219 


-19.3575 








z 5 


-194.829 












Table A23. Same as in Table [All but for the Cousins I c band as 
a function of redshift and the V — I c colour. 





(V-Ic)° 


(y-ic) 1 


(V-Ic) 2 


(y-icf 
















z 1 


-7.92467 


17.6389 


-15.2414 


5.12562 


z 2 


15.7555 


-1.99263 


10.663 


-10.8329 


z 3 


-88.0145 


-42.9575 


46.7401 





z 4 


266.377 


-67.5785 








z 5 


-164.217 












